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Abstract 

The optical spectra of CdSe nanocrystals up to 55 A in diameter are analyzed in a wide range of 
energies from the fine structure of the low-energy excitations to the so-called high-energy transitions. 
We apply a symmetry-based method in two steps. First we take the tight-binding (TB) parameters 
from the bulk sp^s* TB model, extended to include the spin-orbit interaction. The full single- 
particle spectra are obtained from an exact diagonalization by using a group-theoretical treatment. 
The electron-hole interaction is next introduced: Both the Coulomb (direct) and exchange terms 
are considered. The high-energy excitonic transitions are studied by computing the electric dipole 
transition probabilities between single-particle states, while the transition energies are obtained by 
taking into account the Coulomb interaction. The fine structure of the lowest excitonic states is 
analyzed by including the electron-hole exchange interaction and the wurtzite crystal-field terms 
in the exciton Hamiltonian. The latter is diagonalized in the single electron-hole pair excitation 
subspace of progressively increasing size until convergence. The peaks in the theoretical transition 
spectra are then used to deduce the resonant and nonresonant Stokes shifts, which are compared 
with their measured values in photoluminescence experiments. We find that the final results depend 
on the crystal-field term, the relative size of the surface and the degree of saturation of the dangling 
bonds. The results show a satisfactory agreement with the available experimental data. 
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PACS: 71.24+q,71.35.Cc,78.55-m,78.66.Hf 



I. INTRODUCTION 

The optical properties of semiconductor nanocrystals have attracted considerable atten- 
tion in recent years due to their possible applications in quantum dot lasers and other devices 
(see for example Ref. |l|) and, from a more basic point of view, because the nanocrystals are 
the physical realization of small systems where the effect of the low dimensionality should be 
most important. We are interested here in CdSe nanocrystals, which have been intensively 
studied by several groups. iHilli'Si The high-energy excitonic transitions have been investi- 
gated to some extent but most of the experimental works have been devoted to the study 
of the size dependence of the photoluminescence and the fine structure of the low-energy 
exciton states.iflBi In particular, the observation of the "dark exciton" in CdSe quantum 
dots (QD's) is one of the more salient features. 

Several theoretical models have been proposed following two main starting points. First, 
the "particle in the box" point of view is realized in the effective mass approximation (EMA) 
used by several group^'i and the effective-bond-orbital model by Lahel and EinevollJl^ 
From the other point of view, where the QD electron wave function is explicitly constructed 
from the atomic-like orbitals we note the pseudopotentiallHIIll and the semiempirical TlSlli 
methods. 

The TB model that we present here is particularly suitable for describing small nanocrys- 
tals. Comparing to the EMA theories, which are known to be inadequate for the energy gap 
at small sizes, some surface conditions can be varied continuously and their influence in the 
final result evaluated. Also, we analyze the many-body exciton Hamiltonian in a controlled 
way: the electron-hole space is expanded until the convergence in energies is reached. The 
cited EMA theories, in contrast, use a fixed number of electron-hole states. The pseudopo- 
tential calculations published so farffl have been limited to smaller sizes {D < 40 A) than 
ours and some results such as the crystal-field splitting of the valence band seem to be 
controversial: this splitting is expected to be smaller in QD's than in the bulk.lli 

Finally, the previous TB calculation of the exciton fine structure by Leung et al^ is based 
on an unsatisfactory truncation procedure. These authors first use a Lanczos algorithm to 
deduce a few orbital (spin-degenerate) single-particle states below and above the energy gap 
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and then introduce the spin-orbit coupling in the restricted subspace used for diagonalizing 
the exciton Hamiltonian. However, in CdSe nanocrystals, the spin-orbit interaction is almost 
an order of magnitude larger than the average spacing of the TB single-particle levels, and 
drastically modifies the level scheme. Thus their method seems inappropriate beyond the 
immediate vicinity of the exciton ground state. In this paper we propose a different approach 
that allows a rather complete study of the exciton spectrum. It is based on symmetrized 
single-particle states which are first obtained by an exact diagonalization of the full zinc- 
blende TB Hamiltonian including the spin-orbit interaction. The relatively small crystal-field 
term (~ 1/15-th of the spin-orbit one) representing the wurtzite crystal structure is then 
diagonalized in the subspace spanning the exciton states. Not surprisingly, we obtain a 
fine-structure spectrum of the exciton ground state quite different from Ref. 13, especially 
in terms of the relative intensities of the components. Moreover, our study is not limited to 
the lowest-energy excitations: We present results on the nonresonant Stokes shifts and the 
high-energy excitonic transitions. 

This paper is organized as follows. In Section II we give a brief description of the TB 
model, the single-particle states, the exciton Hamiltonian and its diagonalization. In Section 
PI we compare our results with the available experimental data and previous theoretical 
analyses. Finally, in Section |T^ we summarize the main results and present some concluding 
remarks. 



II. TIGHT-BINDING DESCRIPTION OF THE EXCITON STATES 

We consider spherical zincblende crystallites of diameter ranging from 16.63 A (87 atoms) 
to 54.78 A (3109 atoms). In Table | we list some selected crystallite sizes along with their 
number of atoms and number of dangling bonds. The one particle states are obtained 
using the sp^s* nearest-neighbor TB model introduced by Vogl et a/.0 extended here to 
include the spin-orbit interaction. We take the tight-binding parameters for CdSe from Ref. 



r?! , except some minor changes: Es^a = —11-53, Ep^a = 0.53, Eg^c = 1-83, Ep^ = 5.87, 
Vs,s = -3.07, V,,, = 1.8, = 4.23, K.p = 2.17, Vp,s = -5.48, = 7.13, = 6.87, 
Vs*^p = 1.99, Vp^s* = —3.06 eV. The spin-orbit couplings are \se = 0.1434 and Xcd = 0.059 
eV. The dangling bonds are passivated by hydrogen atoms where the H energy level is put 
to Eg^H = —3.3 eV following the scaling prescription given by Kobayashi et al^ We assume 
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that the hopping integrals between the anion or cation and hydrogen follow the Harrison 
scaling rule: Vb-H = {da-c/db-H)^Vac where the db-H are the bond lengths, with b denoting 
cation (6 = c) or anion (b = a). The degree of saturation of the dangling bonds is controlled 
by varying these bond lengths. 

We follow a symmetry-based method developed previously^! to obtain the one-particle 
states. These are computed by diagonalizing the Hamiltonian built in a block-diagonal 
form, in terms of the symmetrized basis corresponding to the Fg, Ti and Fg double valued 
representations of T^. This method allows us to obtain the complete one-particle spectra. 
The eigenstates can be written as, 

0i= E Ci^_fe,„^M^(r-R), (1) 

R,k,m 

where we have omitted b and the symmetry label of the eigenstate for brevity. The ^^(r) 
are the spin-orbit coupled atomic orbitals, following the notation by Koster et a/.0 They are 
constructed from the direct product of the standard basis {s, s* ,Px:,Py,Pz} and the electron 
spin states, all referred to the <001> cubic axis. The s,s* orbitals transform as Fi, the p 
orbitals as F5 and the spin functions as Fg representations of T^. The resulting spin-orbit 
coupled orbitals are thus given hj k = Tq, Fg, F7, and Fg. There are two sets corresponding 
to Cd and Se. The subscript m is the row index and R denotes the atomic site. 

The surface effects can be important in nanocrystals, as pointed out previously0Jli. So, 
we have analyzed here three different cases of dangling bond saturation. First we take 
the bond lengths dcd-H = 1-71 A and dse-u = 1-47 A which correspond to the sum of 
the covalent radii. The analysis of the density of states indicates that there is a relatively 
important presence of surface states near the band edges as in the case of CdTe quantum 
dots.0 The second series is given for dcd-H = 1-21 A and dse-H = 0.97 A. In this case the 
surface contribution has partially disappeared from the states of physical interest. The last 
series of data correspond to dcd-H = 0.81 A and dse-u = 0.57 A; here the contribution 
of the surface states can be found several eV away from the band gap. From these results 
one can assume that the shortening of the db-u can be thought as equivalent to project out 
the surface states from the Hilbert space near the band edge. The physical consequences 
of the dangling bond saturation are also equivalent to the method followed in previous 
TB approaches where the dangling bonds were explicitly removed.lli@ Therefore, we can 
consider the atom-hydrogen bond length as an additional fitting parameter. The study of 
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the excitonic structure has been restricted to the last series because the final results seem to 
be the closest to the measured values. Another surface eff'ect, the surface-to-volume ratio, 
is considered through the number of dangling bonds relative to the total number of atoms 
(cations and anions). In Fig. ^ the distribution of the fraction v = N at / N ctangiing is shown, 
the QD sizes analyzed in this paper are identified by closed squares. Some properties, such 
as the crystal-field splitting of the valence band edge and the photoluminescence Stokes 
shifts seem to be very sensitive to the relative number of dangling bonds. 

In Table we show the size dependence of the one-particle energies and the corresponding 
symmetries of the three highest valence states for the case dcd-H = 0.81, dse-n = 0.57 A. It 
can be seen that the two highest levels are close in energy to each other and well separated 
from the rest of the band for any size. It is interesting to note that this level scheme near 
the 'band edge' fits in with the results of the multiband EMA calculation by Richard et a/.ii 
showing the nearly degenerate quadruplets ISDD3/2 and IPFP3/2. Notice, however, that 



neither of our TB Fg quadruplets is dipole forbidden (see Ref. [191). This energy separation 
of the two fourfold levels from the rest of the band is important. In fact, our numerical 
results for the low-energy exciton states, to be presented later, indicate that the restricted 
e — h subspace spanned by these two highest valence and the lowest conduction levels yields 
almost the same results as larger subspaces. 

The actual crystal structure of the CdSe crystallites is wurtzite.ifl This is taken into 
account as usual by introducing an axial crystal field along a trigonal direction < 111 > 
of the zincblende nanocrystal. This crystal-field splits the atomic p level. In the standard 
sp^s* basis an energy shift has been applied to the pz orbitals.lli In our spin-orbit coupled 
atomic basis the crystal-field Hamiltonian leads to a splitting of the atomic u^^{r) energy 
level into two doublets of C^y symmetry. The corresponding local operator can be written 
as 

i 

+ i i 
-i 1+i 
-i 1-z 

with D = 0.04 eV and J„ the total angular momentum in the < 111 > direction, J„ = 
-^{Jx + Jy + Jz)- This value reproduces the crystal-field splitting of the bulk valence band. 
The matrix elements of Jr, Jy and Jz are obtained in the basis in the standard way.c3 



3 ^ " A' 



D 



(2) 



5 



We take the relation between the total angular momentum basis \JM) and the from 



Ref. 24 



In the basis of the zincblende one-particle states in Eq. (P a matrix element of the total 
crystal-field Hamiltonian is given by, 

{mcM^) = E (C],.rB,™)* C^irB,n(«^'(r - R)|i^c>r(r - R)) (3) 

R,m,,n 

where {u^{r — R)|if"||Mjj*(r — R)) are the matrix elements of the local crystal field Hamil- 
tonian in Eq. (|2|). 

To summarize, the total single-particle Hamiltonian that we use is 

H single = Hq + Hcf, (4) 

where Hcf is given in Eq. (|^) and Hq is the zincblende TB Hamiltonian including the spin- 
orbit interaction. In order to calculate the crystal-field splitting, Ac/, of the highest valence 
band level we diagonalize Hsmgie in the subspace of a sufficiently large number of valence 
states of Hq to reach convergence. The results for several surface conditions are given in 
Table | where it can be seen that the crystal-field splitting presents an irregular behavior 
as a function of the QD size. We find that this behavior is related to the relative number 
of dangling bonds (see Fig. |l]) . In Table | we can see that for a given QD size the splitting 
increases when the hydrogen-atom bond length decreases. This kind of behavior related to 
the surface is also found in the exciton fine structure, even if the non-monotonic jumps are 
not so remarkable as for the crystal-field splitting. 

A complete description of the elementary excitations is made by introducing the Coulomb 
interaction between electron and hole. Let us first define the two-particle states needed 
to describe the excitations above what we call the ground state \g) , describing the fully 
occupied valence band. An electron-hole state can be thought as the pair obtained when 
one electron from the valence band is excited above the gap. The effective hole and electron 
interact through the Coulomb interaction so that we need to describe the excitation in terms 
electron-hole pairs. We call these excited states excitons because they describe the same 
Hilbert state as the bulk excitons when the QD size reach the thermodynamic limit. These 
exciton states can be formally written as 

=^C^^c\v,c), (5) 



where 



\v,c) = ala^lg), (6) 

where the aj(at,) is the creation (annihilation) operator for a conduction (valence) electron 
and \g) is the many-particle ground state. When the Coulomb interaction is introduced the 
matrix elements of the total Hamiltonian can be written in the electron-hole basis after some 
algebra as, 

H^c^v'c' (^c ^v)^vv' ^cc' Jvc^v'c' ~t~ -^vc^v' c' ^ (7) 

with 

Jvcv'c' = (0e(l)0.'(2)|^^^^|0.(2)0c'(l)), (8) 

K^Cv'c' = (0.'(l)0c(2)|^^^^|0.(2)0e'(l)), (9) 

where the (p^ic) are given in Eq. (U), V{r — r') is the bare Coulomb interaction and we 
have explicitly included the dielectric constant e(r — r'). In the preceding expressions it is 
implicitly understood that we consider only excitations involving a single e — h pair. The 
excitations containing two or more pairs and the polarization produced by the surrounding 
ions are taken into account by means of the effective screening in the Coulomb interaction, 
e(r-r'). 

In the Hartree-Fock formulation of the exciton problem adopted here the dielectric con- 
stant is introduced phenomenologically without any distinction between the Coulomb (di- 
rect) and exchange terms. On the other hand, in the many-body formulation in terms of the 
Bethe-Salpeter equationil, the exchange term appears unscreened in the GW approxima- 
tion. Indeed, Rohling and Louie@ have calculated the exciton fine structure in hydrogenated 
Si clusters with an unscreened exchange interaction. However, in the EMA theory of the 
exciton in bulk semiconductors it was early shown0ii that only the short-range part of the 
e-h exchange interaction remains unscreened, but the long-range part is screened. Recently, 
Franceschetti et a/.0 have argued that in a QD the exchange interaction (Eq.(9)) contains 
an important long-range part which needs to be screened. Indeed, they found that within 
their pseudopotential theory, unscreened exchange leads to an excitonic splitting much larger 
than the experimental values. As explained below, we reach a similar conclusion in the TB 
model. There is, of course, no basic controversy over the screening of the Coulomb term. We 
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assume a uniform static dielectric constant, with a size dependence roughly following Ref. 
|30| . As for the screening of the exchange term, we carried out two different calculations: 
(i) completely unscreened, and (ii) unscreened up to the nearest neighbors, but screened 
beyond. A comparison of our results with the experimental Stokes shift data clearly favors 
the second choice. 

The TB approach aims to give an appealing physical description with only a few ad- 
justable parameters. Following this spirit we try to reduce the number of integrals by means 
of some reasonable approximations. Let us take a generic integral in the local spin-coupled 
basis which appears when Eq. (H) and Eq. (||) are expanded in terms of the QD states given 
inEq. ®, 



(n- (r - RiX(r' - K,)\}^^-^\u'^^(r' - Ks)u'^,{v - R4)) (10) 
First, we retain only integrals involving up to two distinct orbitals by setting the sub- 



scripts as either (i) 3=2 and 4=1, or (ii) 3=1 and 4=2. The two choices correspond the 
Coulomb and exchange integrals, respectively. Second, if R2 7^ Ri the Coulomb integral is 
approximated by the monopole-monopole term: l^(|R2 — Ri|). Finally, the on-site Coulomb 
integrals are simplified: Us*s* = Us*p = Us*s = and Uss = Upp = Ugp = U, but U takes 
different values for the cation and the anion. These Coulomb integrals are treated as phe- 



nomenological parameters, but follow qualitatively Ref. ^ for the on-site Coulomb and 
exchange integral values. 



The assumed on-site Coulomb and exchange integrals are given in Table |T|. We have 
checked, however, that the final results are not strongly dependent on these integrals when 
the diagonal values change by 1 eV. Finally, the Coulomb integrals are screened in the 
Coulomb Hamiltonian, they are left unscreened up to the nearest-neighbors (primitive cell 
of the zincblende crystal) in the exchange term and screened otherwise. The on-site screening 
factor is taken to be 0.4 and 0.5 for cations and anions respectively. The nearest-neighbor 
exchange integrals, only important in the exchange Hamiltonian, are assumed to be a tenth 
of the unscreened on-site integrals. 

The exciton Hamiltonian in Eq. (7) is treated by a configuration-interaction-like method 
where we take as many e — h pair states as necessary to get convergence in the final energies 
and energy differences. The size of the exciton space is not a priori fixed. For the high- 
energy transitions, it is adequate to evaluate the average value of the Coulomb term and 
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we take as many as 48 — 50 valence states and 44 — 46 conduction states, depending on 
the QD size. On the other hand, the low-energy fine structure requires a diagonalization 
of the full e — h Hamiltonian including the exchange term and computer time restrictions 
appear; the number of valence and conduction states needed to reach convergence is up to 
18 and 14 respectively. As in previous works0'0 the numerical convergence is faster for 
the energy differences than for the energy levels. The fine structure levels are fixed when 
the energy difference of the first excited level has converged to within 0.1 meV. The highest 
levels analyzed in fine structure present convergence to within 1 meV. 



III. RESULTS 



A. Fine structure 



An early analysis of the band-edge exciton states within the EMA was reported by Efros 
et a/.i; it has been used as a guideline in the TB0 and pseudopotentiallii calculations. In 
the spherical EMA the eightfold multiplet originating from the fourfold highest valence and 
the twofold lowest conduction levels is split by the electron-hole interaction into a fivefold 
passive multiplet and a threefold active triplet. Within our TB model, the symmetry 
analysis leads to three exciton levels, 
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If we call know the exciton Hamiltonian which presents the zincblende symmetry (i.e. 
TB Hamiltonian Hq plus the Coulomb and exchange interaction terms but no crystal field) 
then eigenenergies of these states can be written in terms of the matrix elements of as, 

P _ 1 ttZ _ ttZ I 1 ttZ 

-^3 — 2 3/2,-l/2;3/2,-l/2 3/2,-l/2;-3/2,l/2 2 -3/2,l/2;-3/2,l/2 

Ei = |-f^-l/2 -l/2;-l/2 -1/2 + -^-1/2 -l/2;l/2,l/2 + |-^V2,l/2;l/2,l/2 ' ("'--'-) 
E5 = Es + 2-ff^2,-l/2;-3/2,l/2 

where for the sake of simplicity H^^ m2-m3 m4 stands for (f^i, c^2l-f^^l'^m3? '^rn4}- We compare 
now the analytical result in Eq. ([TTD with the numerical diagonalization of which allows 
us to identify the symmetry of the states. The resulting order for any size is E3 < E4 < E^, 
but E3 and E4 are almost degenerate. 
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The introduction of the crystal-field term reduces the symmetry to C^v, leading to a 
splitting of the previously threefold degenerate levels T^'^ and T^'^. From the Ta to C^y 
compatibility tabled the T'l'^ level is split into T'^'^ + T^'^ and T^^ level into T^'^ + Fg*^. From 
this result we can deduce that the Fg^ singlet is symmetry-forbidden. It is important to 
know if the rest of these first eight levels are truly active. To elucidate this question we 
first project the states obtained from the diagonalization of + Hcf onto the eight states 
obtained from . It is known from group theory that only the T^'^ levels can be optically 
active so we only need the projection of the states obtained from the complete Hamiltonian 
onto these Fg^ active states. 

The resulting projections show that the lowest state, a doublet, is not active and the two 
other doublets and the Fi singlet are symmetry-permitted. However, even this analysis it is 
not enough to conclude that the first doublet is forbidden, it could have other contributions 
from active states of higher energy. We compute then the oscillator strength (see Appendix 
^ for the five different levels. In Table [TVl the detailed results of this analysis are shown for a 
QD of size D = 40.22 A. This allows us to confirm that the first doublet is optically forbbiden 
even when 18 valence and 16 conduction states are included to build the excitonic matrix. 
This result agrees with the EMaS and pseudopotential0 calculations, but not with that of 



Ref. |T3|. We believe this contradiction with the previous TB calculation arises from the 
different procedures followed. We include the spin-orbit interaction from the beginning and 
deduce the full single-particle spectrum. Leung et al. considered this term perturbatively at 
the same level as the Coulomb interaction, applied to a small number of band-edge states 
extracted by using a Lanczos algorithm. 

Let us recall that an analysis of our one-particle spectrum indicates that instead of a 
fourfold degenerate valence level it is more appropriate to consider two fourfold valence 
levels which are very close in energy (see Table ||). Thus a symmetry analysis of the fine 
structure must take into account the first 16 exciton states instead of 8. The analysis of 
this very restricted set of states gives results which are very close to the final values when 
the numerical convergence is reached. The 16 exciton states are structured in general as 
two repeated eightfold multiplets. It is remarkable however that the biggest QD presents a 
special behavior when convergence is reached. In this case the degeneracy pattern of the 16 
states is 2222121121 instead of the 2212122121 for the smaller sizes. 

The resonant Stokes shift is normally compared with the energy difference between the 
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lowest optically active (bright) and the lowest forbidden (dark) exciton states. In Fig. 
H the calculated absorption spectra for several QD's are shown. We identify the first peak 
with the position of the first allowed level. The size dependence of the resulting resonant 
Stokes shift is shown in Fig. | and compared with the available experimental data. There is 
an overall agreement between the theoretical predictions and experimental data. The data 
of the smallest QD's considered here present however a peculiar behavior. The intensity in 
the two smallest QD's in Fig. || shows that the first transition goes to zero and also that an 
additional splitting of the first active state appears and eventually could produce an unclear 
definition of the theoretical Stokes shift. For example, If we associate the mean value of 
the two smallest peaks of the smallest QD studied to the Stokes shift we obtain 28.1 meV. 
On the other hand, one could be tempted to say that the shift measured corresponds to 
the energy difference between the first of the big peaks and the dark level. This feature 
might explain why the experimental values grow so much with decreasing sizes. Another 
interesting feature is the distribution of oscillator strength in EMA. It was shown that in 
spherical and elliptical shape QD's the oscillator strength of the first active state, ±1^, 
relative to that of the 0*^ goes to zero for decreasing QD radius and grows with increasing 
QD size.i We find the same qualitative feature in our model but due to the existence of more 
than three active states a full comparison is not possible. In Fig. |^ there is also another 
interesting aspect concerning the results for two QD's with similar sizes, 30.28 and 31.6 A m 
diameter. In fact, the difference is due to one shell of cadmium atoms. It can be seen that 
the shift increases for the bigger QD, which is opposite to the overall trend of decreasing 
shift with increasing size. We associate this surprising behavior to the relative number of 
dangling bonds in the surface of the QD. This perhaps account for a similar behavior of 
the experimental Stokes shifts. Finally, we have also carried out calculations assuming a 
completely unscreened e-h exchange interaction; the results are shown as open triangles 
connected by dotted lines. The Stokes shift is a little enhanced in small QD's, improving 
the agreement with experimental data. However, in large QD's, the enhancement is too 
large and does not correctly extrapolate toward the bulk exciton splitting. We therefore 
conclude that, at least from a phenomenological point of view, it is necessary to screen the 
exchange interaction beyond the first neighbor, qualitatively agreement with Ref. ^ 

The interpretation of the Stokes shift of the nonresonant photoluminescence seems less 
unambiguous. According to Efros et al, it would correspond to the energy difference between 
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a higher active state and the lowest active one, while Chamarro et al. consider it to be the 
energy difference between a higher active state and the dark ground state. We have adopted 
the latter point of view in presenting the experimental data in Fig. 5. In the TB model we 
analyze the first 16 exciton states. The resulting excited levels, referred to the dark ground 
state, are shown in Figure || along with their relative intensities. The size dependence is 
expressed in terms of the exciton ground state energy as abscissa. One difficulty of the 
comparison between theory and experiment here is that there is no easy way to establish 
an one-to-one relation between our single-QD computed values and the measured points. 
First, the presence of several levels close in energy should be considered. Also, the size 
distribution of the QD's must be considered as was pointed out by Efros et a/.i In Fig. || 
we show again the results for two QD's of similar size as in Fig. 4. Note the significant 
differences between them, in particular, the increase in energy of transition as well as the 
shifting of the absorption intensities to higher energies. Finally, let us mention that, as 
expected, an unscreened e-h exchange leads to larger theoretical values, yielding a better 
agreement with experimental data. However, a fit to the resonant Stokes shift data seems 
more significant. We have therefore retained the screening of the exchange beyond the first 
neighbor. 

B. High-energy transitions 

The high-energy transitions are also analyzed in detail. First the absorption coefficient 
is computed following a simplified procedure (see Appendix ^). Considering the energy 
scale involved, we neglect the exchange term in the interaction Hamiltonian and treat the 
Coulomb term as a first-order perturbation. We identify by inspection the major peaks; 
in Fig. P we show the optical absorption spectra for several QD sizes. We deduce the 
symmetries of the valence and conduction states concerned as well as their relative weights 
when more than one valence-to-conduction transitions is involved. Table shows the results 
for the D = 40.22 A QD. This analysis is done first for the smaller QD, where the peaks 
are well separated and clearly distinguishable. Then the same analysis is performed for 
the bigger QD's, where we try to follow the size evolution of the peaks. The existence of 
level crossing makes it difficult to clearly identify these peaks. For example, the second and 
third peaks for the two smallest QD's are Fg — > F^ and Fg Fg transitions respectively. 
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The two biggest QD's present a reversed order. The evolution of the higher peaks is more 
complicated, due to the mixing of different transitions into each of the peaks. 

The peaks are not immediately identified with the experimental transitions but instead 
in some cases we find it appropriate to group several peaks in the same transition. The 
results are summarized in Fig. ^ 

IV. SUMMARY 

We present a unified picture of the optical properties of spherical CdSe nanocrystals based 
on the sp^s* TB model including the spin-orbit interaction that describes the main features 
of the bulk semiconductor band structure. Our group-theoretical method has allowed us 
to deduce the full one-particle spectra with symmetry-classified eigenstates for crystallite 
sizes up to 55 A in diameter. The dangling bonds were passivated by hydrogen atoms. The 
degree of saturation of the dangling bonds strongly infiuences the electronic properties of 
the QD's due to the presence of surface states near the band edges. The bond lengths from 
the outermost cation or anion to hydrogen were used to study this effect. The final choice of 
the bond lengths removes the surfaces states completely and optimizes the Stokes shift. The 
wurtzite structure of CdSe is treated as usual by introducing a crystal-field term, reducing 
the symmetry from to Cs^. 

With the single-particle states in hand, the exciton states are written in terms of Slater 
determinants, limiting us to the subspace of single electron-hole pair excitations. The 
electron-hole interaction including both the direct (Coulomb) and exchange terms is taken 
into account. In order to derive the fine structure of low-energy excitations, the full exciton 
Hamiltonian is diagonalized in a subspace of progressively increasing size, with as many va- 
lence and conduction states as necessary to reach convergence. The relative electric-dipole 
transition probabilities of the exciton fine-structure components were calculated numerically 
and checked against a symmetry analysis using a restricted subspace. The ground state is 
found to be forbidden in all cases, in agreement with the EMA and pseudopotential cal- 
culations, but in contradiction with the previous TB calculations,^! based on a Lanczos 
algorithm and a perturbative treatment of the spin-orbit interaction. The energy of the 
first allowed peak in the theoretical absorption spectrum yields the resonant Stokes shift of 
the photoluminescence. The theoretical values agree with the measured data except for the 
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very small QD's. The origin of the discrepancy is probably the increasing importance of the 
surface effects in the QD's when their size decreases. Indeed it can be seen that for the small 
QD's {D < 20 A) there are as many or more dangling bonds as atoms. Recently Leung and 
WhaleyB studied the influence of the surface for small QD's and found an increase of the 
Stokes shift after optimized surface relaxation. 

We have also investigated the nonresonant Stokes shift which is associated to allowed 
states lying above the flrst bright level. The theoretical results are systematically smaller 
than the measured values. There are probably several reasons for the discrepancy. First we 
have not included the phonons in our model which seem to be important in the analysis of 
the experimental data (see Ref. | and references therein). Also our results are given for 
each QD individually, we have not taken into account the size distribution of the sample 
which apparently enhances the theoretical shift .ipinally, from a theoretical point of view, 
we find that the exciton fine structure depends rather strongly on the surface conditions: 
both the surface-to-volume ratio and the degree of saturation of the dangling bonds. Also 
the geometry (crystal structure or sample shape) plays an important role, as suggested by 
the increase of the Stokes shifts with increasing crystal field. 

The high-energy excitonic transitions have been also studied. The comparison of the 
theoretical results with the available experimental data shows a reasonable agreement. 
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APPENDIX A: ABSORPTION COEFFICIENT 

The absorption coefficient a is related to €2, a ~ Co'e2,0 where 

^2-^,T.\Mf,?5{E-Ef,) (Al) 
/ 

where iMj^P is squared transition dipole matrix given byEl 

Mf, = {f\[v.H]\g) = -EfJ2C:Ac\r\v) (A2) 
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where Ef is the exciton energy referred to the band gap Eg. It is then easy to check that 
the absorption coefficient can be written, in appropriate units as, 

aiE)r^Y.Ej'\Mfg\'SiE-Efg) (A3) 
/ 

The dipole matrix (e|r|f) is written in terms of the one particles states in Eq. (||) as, 

(A4) 

(«S'(r-Rs')l'5r,K=^(r-R,))], 



where ^r^ = r — . We follow the prescription given in Ref. for the non-zero elements. 
The spin flips are forbidden in optical transitions and we account for that by means of the 
explicit expressions for the on-site and nearest-neighbor dipole matrix elements: 

{u^fy2\5vs\u%^^) = i^die, ^^^^ 



\5rs\u%^) = -iJ\d^e 



+i/2r^si"'+i/2/ — 3 

,6s* |r„ U,8p \ _ ■ f2 



(M_i/2l^r>/3/2) = i^^die, 
and 

{u^,\Srs\u'^,) = 4:d26kk'Smm'ez (A6) 

when u^/ and u'^ belong to nearest-neighbors atoms. In the last case the are only those 
originated by s,px,Py,Pz atomic orbitals. For di and d2 we take the numerical values given 
in Ref. [l^ even though we are only interested in the relative values of the intensity. 

The high-energy transitions have been analyzed by means of a simplified Hamiltonian 
where the exchange term has been neglected. Additionally we take only the diagonal cor- 
rection of the Coulomb term: The expression in Eq. (|A2|) is simplified to (|/) = \vc)), 



Mfg = {f\[r,H]\g) = -(e, - - J,c){c\r\v). (A7) 
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D{k) 


Nat 


^dangling 


0.81 (A) 


1.21 (A) 


1.71 (A) 


16.63 


87 


76 


6.7 


4.3 


3.1 


23.29 


239 


196 


21.0 


20.7 


19.9 


30.28 


525 


276 


16.3 


13.7 


6.1 


31.60 


597 


324 


21.0 


20.6 


17.6 


40.22 


1231 


460 


18.3 


15.3 


8.9 


54.78 


3109 


852 


13.1 


11.5 


7.4 



TABLE I: Some properties of studied QD's are shown. The first column gives the effective diameter, 
the second and third show the number of atoms and danghng bonds, respectively. In the last three 
columns we give the crystal-field splitting (in meV) of the valence band for the different surface 
conditions. 



Diameter 


Vi 


V2 


^^3 


16.63 


-1.034(8) 


-1.055(8) 


-1.163(6) 


23.29 


-0.587(8) 


-0.672(8) 


-0.834(6) 


30.28 


-0.136(8) 


-0.153(8) 


-0.001(6) 


31.60 


-0.380(8) 


-0.413(8) 


-0.578(6) 


40.22 


-0.264(8) 


-0.277(8) 


-0.408(8) 


54.78 


-0.158(8) 


-0.165(8) 


-0.248(8) 



TABLE IL Values of the three highest valence levels (in eV) for different sizes. The symmetry is 
also written in parenthesis. The number in parenthesis 8(6) is short for r8(r6). 
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13(6.5) 


1(0.5) 




1(0.5) 


13(6.5) 



TABLE III: On-site unscreened Coulomb and exchange integrals for cation (anion) in eV. 



Energy (eV) 


Fs fraction 


Oscillator Strength 


2.3076(3) 


0.0 


0.0 


2.3129(3) 


52.47 


1.58 


2.3283(2) 


0.0 


0.0 


2.3307(3) 


43.52 


1.02 


2.3351(3) 


97.29 


1.04 



TABLE IV: Fine structure analysis of the first eight states of the D = 40.22 A quantum QD. 
The first column shows the energies and the corresponding symmetry (in parenthesis). The second 
column shows the relative fraction of F^^^ state (see text) and in the third column we put the 
numerically calculated oscillator strength (in arbitrary units). 



Peak number 


valence-conduction 


Fraction 


1 


8-6 


1 


2 


8-6 


1 


3 


7-6 


1 


1 


8 - 8(8 - 7) 


0.05(0.34) 


5 


7-6(8-6) 


0.99(0.01) 


6 


8-8(6-7) 


0.77(0.17) 


7 


6 - 7(6 - 8) 


0.75(0.25) 



TABLE V: High energy transitions for the D = 40.22 A QD. The symmetry structure of the seven 
first peaks are shown. If more than a transition is concerned the two most important fractions are 
written. 

FIG. 1: Ratio of the number of semiconductor atoms and the number of dangling bonds (open 
circles). The line is the least square linear approximation. The QD's studied in this article are 
shown by closed squares. 
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FIG. 2: The calculated excitonic gap versus QD diameter is compared with the experimental data 
from Ref. |2| (diamonds). The three theoretical curves correspond to different degrees of saturation 
of the dangling bonds, represented by three different sets of cation- and anion-to-hydrogen bond 
lengths (see text). 

FIG. 3: Fine structure optical absorption spectra for several sizes of CdSe nanocrystals. The 
evolution of the first peak, associated to the Stokes shift shows a singular behavior for the smallest 
of the QD's studied here. The size of the two lowest peaks is negligible compared with the next 
peak. This behavior may be related to the sudden increase of the measured Stokes shift in small-size 
QD's (see text). 

FIG. 5: The size dependence of the first excitonic levels (we have omitted the first active level, 
considered before). Note that the abscissa is the exciton ground-state energy, instead of the QD 
size. All the levels are indicated by a cross. The allowed transitions present additionally a closed 
square whose size (area) is proportional to the absorption intensity. The open squares stand for the 
D = 31.60 A case. We also show the measured values of the nonresonant Stokes shifts from Norris 
and Bawendii (open diamonds) and Chamarro et a/.i (open triangles). 

FIG. 6: The size dependence of the absorption coefficient a{E) at high energies. The energy is 
measured from the 'band gap' in each case. 

FIG. 7: High-energy excitonic transition energies are plotted against the the excitonic gap (open 
squares) and compared with the experimental data from Ref. ^ (closed circles). The encircled 
points suggest a probable merger. 

FIG. 4: Resonant Stokes shift obtained from the energy position of the lowest energy peak in Fig. ^ 
(squares), compared with the experimental measurements from Efros et a/.i(crosses) and Chamarro 
et a/.i(closed diamonds). The triangles connected by dotted lines correspond to a completely 
unscreened e-h exchange interaction, while the squares connected by solid lines correspond to a 
finite screening beyond the first neighbors. 
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